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ABSTRACT 

Context. With high-precision radial velocity surveys reaching a sufficiently long time baseline, the domain of long-period planet 
detections has recently opened up. The search for Jupiter-like planets is especially important if we wish to investigate the prevalence 
of solar system analogs, but their detection is complicated by the existence of stellar activity cycles on similar timescales. Radial 
velocity data with sufficiently long-term instrumental precision and robust methods of diagnosing activity are crucial to the detection 
of extrasolar Jupiters. 

Aims. Through our HARPS survey for planets around solar twin stars, we have identified a promising Jupiter twin candidate around 
the star HIPl 1915. We characterize this Keplerian signal and investigate its potential origins in stellar activity. 

Methods. We carry out a Markov chain Monte Carlo (MCMC) analysis of the radial velocity data. To examine the signal’s origin, we 
employ a variety of statistical tests using activity diagnostics such as the Ca II H&K lines and line asymmetry tracers. 

Results. Our analysis indicates that HIPl 1915 hosts a Jupiter-mass planet with a 3800-day orbital period and low eccentricity. Al¬ 
though we cannot definitively rule out an activity cycle interpretation, we find that a planet interpretation is more likely based on a 
joint analysis of RV and activity index data. 

Conclusions. The challenges of long-period radial velocity signals addressed in this paper are critical for the ongoing discovery of 
Jupiter-like exoplanets. If planetary in nature, the signal investigated here represents a very close analog to the solar system in terms 
of both Sun-like host star and Jupiter-like planet. 

Key words, planetary systems - planets and satellites: detection - techniques: radial velocities 


1. Introduction 

As the search for new exoplanets with the radial velocity (RV) 
technique continues, our detection capabilities have increased 
dramatically in two dimensions: smaller signals, such as those 
induced by short-orbit terrestrial-mass planets, can now be found 
thanks to the increased precision of modern RV spectrographs, 
and longer-period signals can also be identified thanks to an ever- 
increasing time baseline of data collection. Both of these detec¬ 
tion capabilities are critical if we wish to find planetary systems 
comparable to our own solar system. 

At the cutting edge of long-period planet detection is the 
search for Jupiter analogs. In terms of signal size, Jupiter in¬ 
duces by far the largest radial velocity shift on the Sun of any 
planet and therefore should be the first planet detected in a twin 
of the solar system. However, the need for more than a decade 
of observational coverage with a consistently precise RV zero 
point complicates the process of discovering such a planet, as 
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does the possibility of stellar magnetic activity cycles inducing 
RV variations on a similar timescale ( |Saar & Donahue|1997[ ). 

In spite of these challenges, a number of Jupiter analog can¬ 
didates have been discovered from long-duration planet surveys 


(see, e.g., Vogt et al.|2000| 

Naef et al.|2003 Wright et al.|2007 

2008 

[Howard et al.[2010 

[Boisse et al. 

2012 

[Marmier et al.| 

2013 

[Wittenmyer et al.|2014 [Feng et al. 

2015 

1 . These surveys 


also demonstrate the prevalence of solar-like magnetic activity 
cycles, which can manifest in the stellar radial velocities ( |Lovis| 
et al.|20lT] Zechmeister et al.|2013|l. These cycles can generally 


be traced by their influence on activity diagnostics such as the Ca 
II H and K lines ( [Baliunas et al.||T995) and asym metries of the 
spectral cross-correlation function ( Dravins|1985[ ). With careful 
attention to these diagnostics, we can minimize the influence of 
activity cycles and recover planet signatures around cycling stars 
( jPumusque et al.|2011| |. 

Over the last four years, we have been conducting an RV 
planet search with the HARPS spectrograph ( [Mayor et al.|2003] l 
targeting a specific subset of 63 stars selected by the proximity 
of their fundamental parameters to those of the Sun (Teff within 
~100 K, log(g) within ~0.1 dex, and [Fe/H] within ~0.1 dex). 
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These solar twin stars offer a unique opportunity to determine 
their compositions at a precision as good as 2% (see, e.g.,|Melen- 


|dez et aL]|2009t [Ramirez et al.||200^ [Bedell et al.||20i4| l. We 
aim to combine this high-precision chemical analysis with the 
m s '-level detection capabilities of HARPS to investigate the 
star-planet connection at an unprecedented level of detail. 

In this paper we present our first detection from the planet 
search program: a Jupiter twin around the solar twin HIPl 1915. 


HIPl 1915 was previously established as a solar twin in Ramfrez 


[et al.|P014[l based on its fundamental properties (Teff = 5760 + 
4 K, log g = 4.46 ± 0.01, and [Fe/H] = -0.059 + 0.004). Its age, 
estimated from isochrone fitting to be 4.0 + 0.6 Gyr, is also in 
keeping with the solar value ( [Ramfrez et al.|20f4) l. HIPl 1915 is 
a bright (my - 8.6) star with relatively low activity (log(R)^^) < 
-4.8), making it an excellent target for high-precision RV moni¬ 
toring. Based on a Keplerian signal in our HARPS RVs, we de¬ 
termine the planet candidate’s orbital parameters. We also inves¬ 
tigate the influence of stellar activity on our RV measurements. 


2. Data 


All of the RV data used in this paper were obtained with the 
HARPS spectrograph on the ESO 3.6 m telescope at La Silla 
Observatory and were processed with the dedicated HARPS 
pipeline. Forty-three data points spanning a time range from Oc¬ 
tober 2009 to January 2015 came from our HARPS large pro¬ 
gram dedicated to solar twin stars (program ID 188.C-0265). An 
additional 13 measurements were taken from the publicly avail¬ 
able archive data and span a time range from October 2003 to 
November 2013 (program IDs 072.C-0488, 089.C-0732, 091.C- 
0034, 092.C-0721, and 183.C-0972). 

For the 43 measurements taken through our program, the 
star was observed in high-accuracy mode with the ThAr lamp as 
a simultaneous reference source. Each observation used a 900- 
second exposure time to minimize potential RV noise from the 
five-minute p-mode oscillations of solar-type stars. To obtain a 
radial velocity measurement, we selected a binary mask based on 
the G2 spectral type for cross-correlation by the pipeline. Owing 
to their generally high signal-to-noise ratios (median S/N over 
the spectral range on the order of 100), the pipeline-produced 


RV error estimates were generally below 1 


We expect 


the actual uncertainties to be somewhat higher owing to stel¬ 
lar noise and potential undiagnosed instrumental effects. To ac¬ 
count for these factors, we added the pipeline-estimated errors 
in quadrature with a baseline error level of 1 m s ', following 
[Dumusque et al.| ( f201 l[ l and in accordance with the scatter in our 
own HARPS survey’s quietest stars. 

The 13 archive measurements have exposure times ranging 
from 120 to 900 seconds. All were taken in high-accuracy mode 
without a source on the simultaneous reference fiber, so that even 
the low S/N measurements are not in danger of spectral contami¬ 
nation by the simultaneous reference source. The data were pro¬ 
cessed by the pipeline and we modified the error estimates as 
described above. An additional 1.1 ms* was added in quadra¬ 
ture to the errors on all exposures under 15 minutes, based on 
the expected p-mode signal for a solar twin ( [Bazot et al.|2011| l. 

We analyzed all of the reduced HARPS spectra to measure 
the Ca II H&K activity indices S hk and log(R/j^) at each obser¬ 
vation epoch. This was done according to the methods ofE ovis[ 
et al. (2011 1 . Errors on Shk were estimated from photon noise. 

Complete RV and S hk data are available in Table along 
with the HARPS pipeline-estimated activity indicators bisector 
inverse span (BIS) and full width at half maximum of the cross¬ 
correlation function (FWHM). 
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3. Analysis and results 


The RMS of the RVs is 6.5 m s ', indicating the presence of 
a signal stronger than the predicted stellar jitter of 2.1 m s~* 
( |Wright|2005| l. Visual inspection and a periodogram analysis of 
the data clearly show a long-period trend in the stellar radial ve¬ 
locities. We ran a grid search across a range of periods from 0 to 
10,000 days and found that the minimum chi-squared value lies 
near 3600 days. 


We inspected the data for signs of stellar activity manifest¬ 
ing in the RVs. A strong correlation between radial velocity and 
the activity index S hk is seen in the measurements from recent 
years, although there is no correlation in the archival data. Sim¬ 
ilar RV-correlated behavior is seen for other activity indicators 
including the bisector inverse slope (BIS), FWHM, Vspan, and 


the biGaussian fit differential (as defined by|Queloz et al.|2001 


[Boisse et al.|201 l[|Figueira et al.|2013] l. Doing a weighted linear 
least-squares fit and removing the trend with any of these activ¬ 
ity indicators from the RVs results in reducing the correlations 
while preserving the long-period planet signal with high signifi¬ 
cance in a periodogram analysis. We additionally attempted the 
prescribed method of Dumusque et al. ( [201 l| l of fitting a long- 
period Keplerian first to an activity diagnostic and then (with 
period and phase fixed) the RVs. This method failed to remove 
the long-period power from the RVs regardless of whether S hk, 
BIS, or FWHM was used as the diagnostic. The failure of these 
methods to successfully “clean” the data of its long-period sig¬ 
nal suggests that the signal originates from a true Doppler shift 
rather than stellar activity, although activity is also present in the 
data. 


3.1. Best model: Keplerian signal + Shk correlation 

In the results presented here, we account for stellar activity by 
fitting a linear trend between RV and S hk simultaneously with 
the Keplerian signal. This way the strength of the correlation be¬ 
tween activity indicator and RV is free to change as the scatter 
induced by the Keplerian signal is removed. We use S hk as the 
primary tracer of activity owing to the unreliability of the BIS di¬ 
agnostic for slow-rotating stars and the demonstrated efficiency 
of Shk in removing activity cycle signals ( [Boisse et al.pOllj 
[Meunier & Lagrange|2013| l. Other potential diagnostics and al¬ 
ternative model choices are explored in Section [T2| 

We ran a Markov chain Monte Carlo (MCMC) analysis to 
fit the radial velocity data to a model consisting of three compo¬ 
nents: a Keplerian signal, a linear correlation term with the S hk 
values, and an offset term. This gives seven free parameters: the 
Keplerian parameters (P, K, e, co + Mq, a> - Mq), where P is 
the orbital period, K is the RV semi-amplitude, e is eccentricity, 
oj is the argument of periastron, and Mq is the mean anomaly 
at reference date fo = 2457000.0; a slope parameter a which 
relates the activity-induced change in RV to the corresponding 
change in S hk', and an offset parameter, C. We include an eighth 
free parameter, cry, representing the jitter that must be added in 
quadrature with the estimated RV error to make the measurement 
uncertainties consistent with the scatter in the residuals. This jit¬ 
ter physically corresponds to unknown Gaussian error that may 
arise from stellar noise or undiagnosed instrumental variations. 
When inputting the data, we subtracted a sigma-clipped median 
value from RV, S hk, BIS, and FWHM for simplicity. This means 
that the offset parameter C is partially set by the difference be¬ 
tween the median RV of the data set and the systemic RV of 
the star, but also includes the RV offset in the linear fit with the 
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Table 1. Best-fit parameters and uncertainties for HIPl 1915b. 


Parameter 


Value 

Uncertainty 

P 

[days] 

3830 

150 

K 

[m s“'] 

12.9 

0.8 

e 


0.10 

0.07 

OJ -1 Mq 

[rad] 

3.0 

1.3 

(jJ - Mq 

[rad] 

2.4 

0.1 

a 

[m s^' (unit5/f;f)^'] 

160 

60 

c 

[m s“'] 

-11.0 

1.3 

o-j 

[m s“'] 

1.8 

0.4 

rrip sin(i) 

[Mjup] 

0.99 

006“^ 

a 

[AU] 

4.8 

0.1 ' 


RMS [m s-i] 29 


Fig. 1. Best-fit orbital solution for HIPl 1915b. Data points are the ob¬ 
served RVs with offset and S hk correlation terms subtracted. One resid¬ 
ual point is located outside of the plot range, as denoted by the arrow. 


median-subtracted S hk values (or any other correlations used in 
the model; see Section|T!2]i. 


All of the fit parameters were given uniform priors except for 
P and K, which were sampled in log-space with a log-uniform 
prior. We ran five parallel chains and used the requirement of 
a Gelman-Rubin test statistic below 1.01 for all parameters to 
avoid non-convergence ( [Gelman & Rubin|1992| l. The best-fit or¬ 
bital solution, plotted in Figure [IJ corresponds to an approxi¬ 
mately Jupiter-mass planet on a 10-year orbit with eccentricity 
consistent with zero (see Table [^. As seen in the MCMC pos¬ 
terior distributions (Figure [^, the chief source of uncertainty in 
the period measurement is in its degeneracy with eccentricity 
and the associated argument of periapse, parameters that are in¬ 
evitably poorly constrained when only one orbit has been ob¬ 
served. The best-fit jitter term is consistent with the pred icted 


stellar jitter of 2.1 m s ^ from the formula of Wright (2005 i. 


After subtracting the best-fit planet signal, we searched for 
additional signals in the residuals with a generalized Lomb- 
Scargle periodogram analysis ( [Zechmeister & Kurster|200^ . No 
significant periodicities were identified in the residuals to the 
planet fit or in the residuals to the planet and the fitted S hk trend. 
We place constraints on the presence of additional planets in the 
system using a test based on the bootstrapping and injection of 
a planet signal into the residuals. For this test, we use the resid¬ 
uals to the planet fit only, and inflate the error bars (with jitter 
included) such that a fiat-line fit to the residuals gives a reduced 
chi-squared of one. We then run a set of trials in which the RV 
data are resampled with replacement at the times of observa¬ 
tion, a planet signal with period and semi-amplitude fixed and 
the other orbital parameters randomized is injected, and the re¬ 
duced chi-squared of a flat-line fit to these simulated data is mea¬ 
sured. If the chi-squared statistic has increased with three-sigma 
significance in 99 out of 100 trials at a given period and velocity 
semi-amplitude, we consider such a planet signal to be excluded 
based on its inconsistency with the scatter in the actual RVs. We 
run this trial until upper limits on planet presence are found for 
all of the periods in a log-sampled grid. The results rule out the 
presence of gas giants within a 1000-day orbital period, leaving 
open the possibility of inner terrestrial planets (Figure]^. 


a (A.U.) 

0.1 1.0 



Fig. 2. Exclusion limits on potential planets in the RV data after subtrac¬ 
tion of the Keplerian signal. Any signal in the parameter space located 
above the lines plotted would have induced a larger scatter in the RV 
residuals with three-sigma significance if present. Potential trends with 
activity indicators were not removed for this test to avoid underestimat¬ 
ing the RV scatter. 


3.2. Alternative model choices 


It is not obvious from theory that the Keplerian + S hk correla¬ 
tion model used above should be the best model choice. Other 
activity indicators exist and have been used with success to re¬ 
move activity trends from RVs in the past. Chief among these are 
BIS and FWHM, both characteristics derived from the RV cross¬ 
correlation function (CCF) which measure the average spectral 
line asymmetry and/or broadening at the epoch of the observa¬ 
tion (see, e.g., |Queloz et al.|2001|[2009] l. Either BIS or FWHM 
could potentially be used in place of S hk- Also, because these in¬ 
dicators arise from distortions to the rotationally broadened stel¬ 
lar absorption lines while S hk probes chromospheric emission, 
we could combine S hk with a CCF-derived measurement to si¬ 
multaneously track two separate physical indicators of activity. 

We considered several alternative models before settling on 
the Keplerian + S hk model as the best choice. A summary of the 
models tried is given in TableIn all cases, the fitting was done 
with an MCMC as described above. The inclusion of a Keplerian 
signal is strongly justified, with the Keplerian H- 5 //a- fit having 
a false alarm probability on the order of 10 compared to the 
best fitting activity model (S hk + FWHM). While the lowest 
value is obtained for the fit to a Keplerian, S hk, and FWHM, the 
addition of the FWHM term is not fully justified by the reduc- 
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Table 2. Models considered. 


Name 

Degrees of Freedom 

of Fit 

Keplerian only 

50 

162.54 

Keplerian + S hk 

49 

136.00 

Keplerian H- FWHM 

49 

145.32 

Keplerian + BIS 

49 

133.96 

Keplerian + S hk + FWHM 

48 

129.50 

S hk only 

54 

897.78 

BIS only 

54 

704.27 

FWHM only 

54 

650.97 

Shk + BIS 

53 

719.28 

Shk + FWHM 

53 

639.76 


tion in scatter around the fit. An F-test gives a formal false alarm 
probability of 13% for the Kepleiian + Shk + FWHM model 
compared to the simpler Keplerian + Shk model. The Keplerian 
+ fit has a false alarm probability of 0.3% when compared 
to the Keplerian only model. Using BIS in place of S hk techni¬ 
cally gives a slightly better fit, but we chose to use S hk in the 
main analysis because of BIS’s previously noted unreliability for 
slow rotators. Regardless of the choice of activity indicator(s), all 
Keplerian fit parameters remain within the uncertainty intervals 
quoted in Table [T] 


4. Planet or activity cycle? 


The question of the role of stellar activity role in the observed 
radial velocities is of special importance given the proximity of 
the signal period to the 11-year solar magnetic cycle. The sig¬ 
nal size is also consistent with the theor etically predicted 10- 
11 m s“' signal for a solar activity cycle (Meunier & Lagrange 
|2013| l. However, an activity cycle should manifest not only in the 
observed RVs, but also in other diagnostics. We expect a pos¬ 
itive correlation in RV with the activity index S hk as well as 
with other diagnostics such as BIS and the FWHM of the cross¬ 
correlation function ( |Lovis et al.|2011[|Zechmeister et al.|2013| l. 
While such correlations are seen, they appear insufficient to fully 
explain the observed RVs. 

For at least two other solar twins in our sample, similar 
Jupiter-like RV variations are seen; however, in these cases sub¬ 
tracting a linear fit with S hk entirely removes the long-period 
power from a periodogram and reduces the RMS of RVs to a 
level consistent with the expected stellar jitter. In contrast, the 
RV signal inferred from the variation of S hk or any other activ¬ 
ity tracer for HIPl 1915 is much smaller in size than the total RV 
signal, so that significant power remains at long periods in the 
periodogram. 

A large sample of stars with activity cycles have been stud¬ 
ied with HARPS by |Lovis et ah] ( |201 1[ ), who derived empirical 
predictions for the expected correlation between RV and other 
activity indictors with as a function of stellar temperature 
and metallicity. The correlations seen for HIPl 1915 are gener¬ 
ally lower than the predicted values: Crv = 30 + 1 m s ' (unit 
Cfwhm = 58 + 2 m s"* (unit and Cbis = 22.4 

+ 0.7 m s“' (unit R'hk) ^ from their Equations 9-12, while the C 
values recovered from our data are 11, 45, and 16 respectively. 
This disagreement suggests that the strength of activity seen in 
the data is insufficient to explain the RV signal. 

An alternative interpretation is that we are observing activity 
overlaid with a planet signal in the radial velocities. In this view, 
the correlation with activity indicator(s) in the RV residuals is 
independent of the Keplerian signal. One piece of evidence in 
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favor of this interpretation comes from inspecting the correlation 
plots before and after removal of the 3800-day Keplerian signal. 
If this signal is truly unrelated to the stellar activity, it should 
add random scatter about the correlation. Indeed, subtracting the 
Keplerian signal from the data reduces the scatter about the best- 
fit trend between RV and Shk from 7.7 to 3.0 m s ' and has a 
comparable effect on the correlations with other activity indices 
(Figure [^. 

Inspecting the behavior of S hk with time sheds some light on 
the evolution of HIPl 1915’s activity level. As depicted in Figure 
1^ the stellar activity level exhibited an incoherent scatter un¬ 
til mid-2013, at which point it increased from the equivalent of 
log(Rj^^)—4.87 to —4.82 by mid-2014. A concurrent rise in 
BIS and FWHM is also seen at this epoch. This level of vari¬ 
ability is consistent with the behavior of the Sun, which has a 
mean log(Rj^^) of -4.94 and an amplitude of variation around 
0.1 dex over a full activity cycle ( [Hall et ar]|2007[ ). Given the 
recent increase in Ca emission, we could then expect a corre¬ 
sponding increase in short-duration activity phenomena such as 
spots and plages to manifest in the more recent RVs. This would 
add power to the correlations with line asymmetry indicators. 

Ideally, we would be able to remove activity phenomena such 
as spots and plages by averaging our measurements over the 
timescale of the stellar rotation period, on the order of a month. 
Unfortunately the sampling of the data is too uneven to aver¬ 
age effectively over the rotation period. We average instead over 
an entire observing season, which is a regime in which the time 
coverage should be distributed evenly enough to remove activity 
signals at periods shorter than a cycle. The resulting data from 
this seasonal binning trace the 3800-day signal well, but show 
no statistically significant correlations with activity indicators. 

In short, we have considered a variety of activity diagnostics 
and methods of cycle removal and the 3800-day signal remains 
robustly present in the data. The coincidence in phase between 
the rise in S hk and the 3800-day signal maximum, however, 
calls for further monitoring of the star before we can definitively 
rule out the possibility of an activity cycle inducing the RV sig¬ 
nal. 

5. Summary 

We have detected a Keplerian signal corresponding to a gas gi¬ 
ant planet with a 10-year orbital period, consistent with a Jupiter 
twin, around the solar twin HIPl 1915. While we cannot conclu¬ 
sively exclude an activity cycle interpretation, commonly em¬ 
ployed activity tracers fail to account for the signal. If the 3800- 
day signal is entirely induced by an activity cycle, the failure of 
various techniques to remove it bodes poorly for the prospect of 
correcting for such cycles in future long-period planet searches. 

If this signal is truly planetary in origin, the HIPl 1915 sys¬ 
tem is a close analog to the solar system. Our analysis shows the 
planet HIPl 1915b to be a close match to Jupiter both in mass 
and orbital period, and its host star is extraordinarily similar to 
the Sun. Beyond having fundamental physical properties close 
to the Sun, initial spectroscopic analysis suggests that HIPl 1915 
is also a solar twin in the sense that its detailed abundance pat¬ 
tern matches the solar pattern (Melendez et al. 2015, in prep). 
The presence of a Jupiter twin and a solar-like composition both 
make HIPl 1915 an excellent prospect for future terrestrial planet 
searches. 
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Fig. 3. Correlations for the activity tracers 5 hk, BIS, and FWHM with radial velocity. The upper panels shows the observed RVs. The lower panels 
show the residuals to the 3800-day signal fit (as presented in Section[3T]l. Dashed lines are linear least-squares best fits to the data. 
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Fig. 4. Evolution of the activity indices Shk, BIS, and FWHM with time. Errors on BIS and FWHM are approximated as the errors on the RVs. 
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Fig. 5. Constraints on the fit to the radial velocities from the MCMC analysis. Parameters are defined in Section Plot made using code from 
|Foreman-Mackey et al.|([^14|>. 
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Table 3. HARPS measured radial velocities and activity indices for 
HIPl 1915. 


Julian Date 

Exp. Time (s) 

RV (km s-‘) 

ctrv (km s 

Shk 

0 ' s ™ 

BIS (km s ') 

FWHM(kms-0 

2452941.7219 

116.5 

14.5766 

0.0030 

0.1856 

0.0020 

-0.0200 

6.7670 

2453271.8214 

182.9 

14.5622 

0.0038 

0.1987 

0.0027 

-0.0358 

6.8019 

2453312.6877 

900.0 

14.5722 

0.0017 

0.1691 

0.0013 

-0.0227 

6.7696 

2453987.7974 

121.3 

14.5610 

0.0021 

0.2073 

0.0017 

-0.0260 

6.7872 

2454438.6089 

134.1 

14.5571 

0.0019 

0.1829 

0.0015 

-0.0364 

6.7765 

2454752.7519 

300.0 

14.5529 

0.0015 

0.1915 

0.0013 

-0.0312 

6.7720 

2454764.7361 

900.0 

14.5507 

0.0016 

0.1853 

0.0013 

-0.0377 

6.7782 

2455217.5801 

121.9 

14.5537 

0.0023 

0.1983 

0.0018 

-0.0355 

6.7863 

2455543.6010 

500.0 

14.5633 

0.0015 

0.1738 

0.0012 

-0.0308 

6.7929 

2455846.8020 

900.0 

14.5692 

0.0011 

0.1868 

0.0010 

-0.0327 

6.7871 

2455850.7947 

900.0 

14.5698 

0.0009 

0.1854 

0.0009 

-0.0328 

6.7891 

2455851.7841 

900.0 

14.5697 

0.0008 

0.1812 

0.0008 

-0.0285 

6.7817 

2455852.7779 

900.0 

14.5697 

0.0010 

0.1850 

0.0009 

-0.0328 

6.7861 

2456152.8532 

243.8 

14.5660 

0.0020 

0.1694 

0.0016 

-0.0394 

6.7951 

2456164.9273 

900.0 

14.5716 

0.0009 

0.1814 

0.0008 

-0.0317 

6.7822 

2456165.9242 

900.0 

14.5721 

0.0009 

0.1810 

0.0008 

-0.0311 

6.7800 

2456299.6151 

900.0 

14.5749 

0.0012 

0.1870 

0.0011 

-0.0355 

6.7922 

2456300.5833 

900.0 

14.5750 

0.0006 

0.1815 

0.0006 

-0.0298 

6.7854 

2456486.9106 

900.0 

14.5774 

0.0007 

0.1834 

0.0007 

-0.0309 

6.7898 

2456487.8968 

900.0 

14.5778 

0.0007 

0.1850 

0.0007 

-0.0328 

6.7917 

2456488.9425 

900.0 

14.5749 

0.0008 

0.1800 

0.0007 

-0.0308 

6.7857 

2456489.9048 

900.0 

14.5773 

0.0008 

0.1840 

0.0008 

-0.0266 

6.7862 

2456490.8972 

900.0 

14.5753 

0.0007 

0.1789 

0.0007 

-0.0334 

6.7856 

2456502.9393 

900.0 

14.5796 

0.0023 

0.1672 

0.0017 

-0.0371 

6.7938 

2456557.9085 

900.0 

14.5740 

0.0020 

0.1825 

0.0016 

-0.0326 

6.8084 

2456558.8337 

900.0 

14.5793 

0.0007 

0.1892 

0.0006 

-0.0303 

6.7940 

2456559.8500 

900.0 

14.5819 

0.0007 

0.1875 

0.0007 

-0.0280 

6.7945 

2456560.8411 

900.0 

14.5831 

0.0007 

0.1884 

0.0006 

-0.0302 

6.7946 

2456589.8019 

400.0 

14.5790 

0.0010 

0.1902 

0.0010 

-0.0308 

6.8011 

2456615.6781 

300.0 

14.5798 

0.0024 

0.1969 

0.0019 

-0.0259 

6.7981 

2456637.6474 

900.0 

14.5789 

0.0007 

0.1884 

0.0007 

-0.0256 

6.7904 

2456638.6662 

900.0 

14.5815 

0.0008 

0.1913 

0.0008 

-0.0279 

6.7949 

2456639.6356 

900.0 

14.5802 

0.0008 

0.1866 

0.0007 

-0.0293 

6.7955 

2456640.6347 

900.0 

14.5783 

0.0008 

0.1904 

0.0007 

-0.0258 

6.7975 

2456641.6456 

900.0 

14.5786 

0.0007 

0.1862 

0.0007 

-0.0284 

6.7921 

2456643.6517 

900.0 

14.5778 

0.0010 

0.1877 

0.0008 

-0.0267 

6.7923 

2456708.5010 

900.0 

14.5764 

0.0009 

0.1915 

0.0009 

-0.0255 

6.8008 

2456709.5310 

900.0 

14.5795 

0.0008 

0.1917 

0.0008 

-0.0312 

6.7933 

2456710.5234 

900.0 

14.5798 

0.0009 

0.1951 

0.0009 

-0.0242 

6.7964 

2456711.5226 

900.0 

14.5796 

0.0007 

0.1916 

0.0008 

-0.0253 

6.7965 

2456851.9220 

900.0 

14.5807 

0.0009 

0.1993 

0.0007 

-0.0259 

6.8022 

2456852.8915 

900.0 

14.5823 

0.0012 

0.1949 

0.0011 

-0.0223 

6.7973 

2456855.8977 

900.0 

14.5760 

0.0011 

0.1981 

0.0011 

-0.0312 

6.7999 

2456856.8881 

900.0 

14.5783 

0.0010 

0.1923 

0.0009 

-0.0283 

6.8029 

2456904.8198 

900.0 

14.5825 

0.0010 

0.1936 

0.0009 

-0.0266 

6.8065 

2456906.7954 

900.0 

14.5812 

0.0010 

0.1969 

0.0010 

-0.0262 

6.8068 

2456907.8079 

900.0 

14.5785 

0.0009 

0.1934 

0.0009 

-0.0264 

6.8062 

2456960.7641 

900.0 

14.5750 

0.0007 

0.1876 

0.0007 

-0.0291 

6.7933 

2456961.7200 

900.0 

14.5726 

0.0007 

0.1873 

0.0006 

-0.0306 

6.7951 

2456963.7068 

900.0 

14.5722 

0.0007 

0.1891 

0.0006 

-0.0310 

6.7938 

2456964.7840 

900.0 

14.5737 

0.0007 

0.1887 

0.0007 

-0.0271 

6.7961 

2456965.6863 

900.0 

14.5743 

0.0006 

0.1909 

0.0006 

-0.0273 

6.7971 

2456966.7461 

900.0 

14.5757 

0.0009 

0.1939 

0.0008 

-0.0244 

6.7943 

2457025.5721 

900.0 

14.5755 

0.0009 

0.1897 

0.0008 

-0.0269 

6.8001 

2457026.6043 

900.0 

14.5745 

0.0007 

0.1896 

0.0006 

-0.0266 

6.7993 

2457027.5960 

900.0 

14.5727 

0.0008 

0.1921 

0.0007 

-0.0261 

6.7989 
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